Suppression of low-energy Andreev states by a supercurrent in YBa 2 Cu 3 07 s 



in 
o 
o 



oo 

(N 
O 

o 



^3 

G 

O 

o 



OS 

SO 
in 
o 
in 
o 

-I— > 



i 

c 

o 
o 



X 



J. Ngai, P. Morales, and J.Y.T. Wei 

Department of Physics, University of Toronto, 60 St. George St., Toronto, ON M5S1A1, Canada 

(Dated: February 2, 2008) 

We report a coherence-length scale phenomenon related to how the high-T c order parameter (OP) 
evolves under a directly- applied supercurrent. Scanning tunneling spectroscopy was performed on 
current-carrying YBa2Cu307_s thin-film strips at 4.2K. At current levels well below the theoretical 
depairing limit, the low-energy Andreev states are suppressed by the supercurrent, while the gap-like 
structures remain unchanged. We rule out the likelihood of various extrinsic effects, and propose 
instead a model based on phase fluctuations in the rf-wave BTK formalism to explain the suppression. 
Our results suggest that a supercurrent could weaken the local phase coherence while preserving 
the pairing amplitude. Other possible scenarios which may cause the observed phenomenon are also 
discussed. 

PACS numbers: 74.72.-h, 74.40.+k, 74.45.+C, 74.50.+r 



I. INTRODUCTION 

High critical-temperature (T c ) cuprate superconduc- 
tors are distinguished from conventional s-wave super- 
conductors by a predominantly <i-wave order parame- 
ter (OP)i£. The rf-wave OP symmetry has led to a 
wealth of unconventional experimental phenomena. One 
such phenomenon is the formation of low-energy An- 
dreev states bound on non principal-axis surfaces of 
cuprate superconductors^*^. When probed by quasi- 
particle tunneling, these Andreev states are manifested 
as zero-bias peaks (ZBCP) in the conductance spec- 
tra of high-impedance tunnel junctions. These Andreev 
states form when time-reversed quasiparticles interfere 
constructively via the phase sign change of the d-wave 
OP about its nodal (k x ~±k y ) axes&i For s-wave su- 
perconductors, tunneling involves quasiparticle transmis- 
sion, and is thus only sensitive to the amplitude of 
the OP. Tunneling into c?-wave superconductors involves 
both quasiparticle transmission and Andreev reflection, 
and is thus sensitive to both the phase and amplitude of 
the pair wavefunction&^i. 

Recent interest has been focused on how the d- 
wave Andreev bound states, quasiparticle density of 
states, and pairing symmetry itself could evolve un- 
der an applied supercurrenljiSiSiifi. In this paper, 
we present a scanning tunneling spectroscopy (STS) 
study of YBa2Cu307_5 (YBCO) thin-film strips carry- 
ing directly-applied supercurrents at 4.2K. Our study 
uses the supercurrent to perturb the superconducting 
OP electrodynamically, rather than thermodynamically 
as was previously done in magnetic field studiesii*^. 
High- T c cuprates have characteristically short coherence 
lengths, which can be ~nm in the basal plane. By 
virtue of its nanoscale junction size, STS can thus provide 
coherence- length scale information about the high- T c OP 
when driven by an applied supercurrent. In our experi- 
ment, we observed systematic suppression of the phase- 
sensitive spectral features by the supercurrent, while the 
amplitude-dependent features remain largely unchanged. 
We discuss and rule out the likelihood of various ex- 



trinsic effects, and propose instead a model based on 
current-driven OP phase fluctuations to describe the phe- 
nomenon observed. 



II. EXPERIMENTAL SETUP AND RESULTS 

Our experimental setup is illustrated in the left in- 
set of Figure ^ showing a Pt-Ir tip positioned over a 
YBCO thin-film strip in an STS geometry, with a Au 
contact providing the bias voltage V between tip and 
sample. STS is performed in the usual manner, by mea- 
suring the tunneling current I t with the scanning and 
feedback temporarily suspended, but while a current I s 
is simultaneously applied through the superconducting 
film at 4.2K. In this STS geometry, the high junction 
impedance (~0.1G£!) ensures that I t (~lnA) is decou- 
pled from I s (~100mA). The I s is supplied from a float- 
ing source below, via Au contacts on the ends of the 
film strip. A special circuitry is used to synchronize I s 
with the STS feedback and data acquisition 13 . Short 
duty-cycle I s pulses «200/is in width are used to pre- 
vent sample Joule heating, which could occur at high 
current levels. DC measurements are also made at lower 
current levels for comparison, in order to rule out tran- 
sient electronic effects. For our experiment, this versatile 
STS technique enabled very low-noise (<5pA) and highly 
reproducible tunneling I-V measurements at nanometer 
length scales. 

The YBCO thin-film strips used in our experiment 
were grown by pulsed laser-ablated deposition on SrTi03 
(STO) substrates. Both {110} and {001} oriented films 
were made, the former using a templating technique. The 
films were typically 50nm thick and in a 1mm x 3mm 
strip-line geometry, aligned along <100> and [110] re- 
spectively for the {001} and {110} films. The film strips 
showed well-defined T c 's («87K) and sizable transport 
critical-current densities at 4.2K ( J c «5x 10 6 A/cm 2 , de- 
fined by the appearance of 1/iV/mm along the strip). 
Epitaxiality of the films was confirmed by x-ray diffrac- 
tion and rocking-curve analysis. Film surface quality was 
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FIG. 1: (Color online) Normalized dl/dV spectra obtained on 
a {110} YBCO film strip carrying a supercurrent / s applied 
along [110] at 4.2K. Right inset displays the unnormalized 
data. Left inset displays our experimental setup showing an 
STM tip over a YBCO thin-film strip grown on an STO sub- 
strate with the orientation as indicated. The tip is biased 
relative to the sample by a gold pad on the sample surface. 
Spectroscopy measurements are made while a simultaneous 
current is passed through the film strip. 



determined directly from topography measurements us- 
ing our STS setup. Shown in the inset of Figure is 
a topographic image taken on a {001} film, indicating 
smooth terraces with large c-axis faces ~30nm in size 
and a6-plane edges ~3nm in height. 

Figure ^ shows tunneling conductance data taken at 
4.2K on a {110} film strip at various current levels ap- 
plied along [110]. The dl/dV spectra, which were nu- 
merically differentiated from the I-V data, are plotted 
in the right inset, and the normalized dl/dV spectra are 
plotted in the main panel. These spectra show a well- 
developed ZBCP structure which is suppressed in height 
but only slightly changed in width by the applied current. 
Such ZBCP structures have been commonly observed in 
YBCO, and could be identified as the low-energy An- 
dreev states based on a predominant rf-wave OB2^&. 
Two features of the ZBCP spectral evolution versus cur- 
rent are noteworthy. First, the ZBCP shows no splitting, 
indicating time-reversal invariance of the d-wave OP up 
to 160mA (the maximum applied current) in our near- 
optimally doped YBCO samples at 4.2K. Second, the 
area under dl/dV is not conserved, indicating negligible 
thermal broadening, with a clear loss of Andreev states. 
The non-splitting of the ZBCP we observed is quantita- 
tively consistent with a recent theoretical calculation? , 
which shows no Doppler shiftii*i£ of the Andreev bound 
states in this geometry until a higher supercurrent den- 
sity. However, the clear loss of Andreev states we ob- 
served was not anticipated by this mean-field calculation. 



III. ANALYSIS 

We have considered various extrinsic effects in analyz- 
ing the phenomenon we observed. One possible extrin- 
sic effect is the Joule heating from the sizable current 
applied. However, our pulsed current technique, along 
with the absence of significant ZBCP broadening, largely 
rule out this possibility. A second possible extrinsic ef- 
fect involves current-generated vortices that could intro- 
duce Aharonov-Bohm phases in the Andreev interference 
process^. We note that our current levels are well be- 
low the theoretical depairing limrfci, and also below any 
sharp voltage onsets associated with vortex depinning, 
although our present experimental resolution cannot yet 
rule out the presence of such vortices. Nevertheless, the 
theoretical calculation in Ref.14 has predicted a ZBCP 
broadening which we do not observe, thus rendering this 
vortex scenario unlikely. A third possible extrinsic effect 
involves excess low-energy quasiparticles due to the ap- 
plied nodal current^. In this scenario, the excess quasi- 
particles could effectively raise the quasiparticle tempera- 
ture, again creating spectral broadening which we do not 
observe. In summary, all three extrinsic effects described 
above would be expected to yield a conservation of the 
low-energy Andreev states, and are therefore unlikely to 
explain our observations. 

Having ruled out the likelihood of these various extrin- 
sic scenarios, we propose instead an explanation based 
on current-driven OP phase fluctuations to describe the 
observed phenomenon. High-temperature superconduct- 
ing cuprates are distinguished by their low carrier den- 
sities and short coherence lengths^. These proper- 
ties imply an inherently low superfluid stiffness, mak- 
ing the order parameter susceptible to fluctuations in 
its phaseit Motivated by this possibility, we incor- 
porate phase fluctuations into the generalized Blonder- 
Tinkham-Klapwijk (BTK) 18 formalism, which has been 
theoretically established 4 for rf-wave superconductors, 
and experimentally verified for YBCO^. The generalized 

BTK expression for tunneling conductance G ns is given 
by^4& 



where a and b are the Andreev-reflection and normal- 
reflection coefficients, f is the Fermi-Dirac function, G nn 
is the normal-state junction conductance, E the quasi- 
particle energy, and 9 the polar angle in fc-space. To 
model the effects of OP phase fluctuations, we modify 
the BTK kernel from Ref.4 by adding a phase factor e %v 
to the phase-interference term in its denominator (see 
Eqn|21below) then integrating it over the domain of phase 
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FIG. 2: (Color online) Phase fluctuation model showing the 
evolution of the ZBCP on a {110} junction as a function of 
the phase smearing parameter 7. Left inset illustrates the 
Andreev interference process based on the d-wave phase sign 
change. Right inset shows the ZBCP height and spectral area 
as a function of 7. 
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Note that the amplitude factors T± — (i£/|A±|) — 
\J (_E/|A±|) 2 — 1 contain a d-wave gap function with Ao 
as gap maximum. The phase factors e l ^ ± = A±/|A±| 
represent the sign of the pair potential A± = A(0±) 
which is experienced by an Andreev-reflected quasielec- 
tron (or quasihole) propagating at an angle 8 + (or 6- = 
7r — 6+) relative to the junction normal (see left inset in 
Fig. 3). Through the extra phase introduced by e lv , the 
weighted integral over ip essentially smears out the rela- 
tive phase between consecutive Andreev-reflected quasi- 
particles. This phase smearing would cause incoming 
and outgoing quasiparticles to interfere less construc- 
tively, thereby leading to a loss of the low-energy Andreev 
states. Our phenomenological approach can be compared 
with the more rigorous treatment of Ref.19 which uses a 
2D X-Y model based on vortex fluctuations^,. Our model 
makes no assumptions about which physical mechanism 
is actually driving the phase fluctuations, while also go- 
ing beyond the low-impedance (Z<^1) and principle-axis 
({100}) junction cases^. 

Our spectral simulations for a {110} junction are plot- 
ted in Figure |2 as normalized conductance G ns /G nn . 
A large barrier parameter Z—1.5 was used to represent 
our high-impedance junctions, along with the experimen- 
tal temperature T=4.2K and a typical gap maximum 
Ao=15meV. A clear trend of ZBCP suppression is seen 



FIG. 3: (Color online) Normalized dl/dV spectra obtained on 
a {001} YBCO thin-film strip carrying a supercurrent I s ap- 
plied predominantly along <100> at 4.2K. Right inset shows 
the evolution of the ZBCP height and spectral area as a func- 
tion of the applied current. 



versus the phase-smearing parameter 7, showing a drop 
in the peak height and a loss of spectral area. Detailed 
evolutions of height and area versus 7 are plotted in the 
right inset of Fig|5J Comparing our data in Fig^with the 
simulations in Figure|21 it is clear that the ZBCP suppres- 
sion could be qualitatively attributed to the dephasing 
effects of current-driven OP phase fluctuations. Physi- 
cally speaking, consecutive Andreev-reflected quasiparti- 
cles which are bound within a coherence length £ from 
the specular {110} surface (see left inset of Fig. 2), would 
experience an effective smearing in their relative phase 
as a result of the phase fluctuations, thus interfering less 
constructively. This phase smearing leads directly to a 
loss in the low-energy Andreev states, resulting in the 
ZBCP suppression seen in the data. 



IV. DISCUSSION 

It is important to consider possible physical mech- 
anisms by which the supercurrent is driving the local 
phase fluctuations. One scenario involves the possibil- 
ity of spatial doping inhomogeneities in the supercon- 
ducting cuprates 21 . In the presence of current flow, lo- 
cal areas of strong superconductivity could become de- 
coupled from each other, thus permitting local phase 
fluctuations. Another scenario may involve the prolif- 
eration of thermally-generated vortex-antivortex pairs in 
the Kosterlitz-Thouless pictur e) 22 ! 23 , but with current as 
an additional driving forc o 24 i 25 . In principle, the static 
phase gradient carried by the supercurrent could also 
couple dynamically with other types of phase-ordered 
states, such as spin/charge-density wave a 26 ' 27 , staggered 
flux- like states 2 ^ 2 ^ 3 ^ or collective mode a 31 i 32 , to produce 
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FIG. 4: (Color online) Normalized dl/dV spectra obtained 
on a {001} YBCO thin-film strip carrying a supercurrent I s 
applied predominantly along <100> at 4.2K. The inset shows 
film topography. 

enhanced phase fluctuations. In general, any of the above 
scenarios may lead to local regions containing a weak- 
ened OP, with respect to phase fluctuations. In this con- 
nection, it is worth noting that such weakened regions 
could also be nucleation centers for phase-slip phenom- 
ena, which were seen in recent transport studies on su- 
perconducting YBCO microstripsiS^i. 

It should also be remarked that our discussion of low- 
energy Andreev states is not limited to only specular 
cases. These Andreev states could also form under non- 
specular conditions in a d-wave superconductor, i.e. near 
grain boundaries, extended non-magnetic impurities, or 
surface interstitial defecta 35 i 36 i 37 i 38 i 39 i 40 i 41 . Like the spec- 
ular Andreev bound states, these non-specular resonance 
states are essentially also based on rf-wave Andreev in- 
terference, and should therefore also be conserved under 
an applied current, as was shown in a recent theoretical 
calculation 42 . However, another calculation for this non- 
specular case, done as a function of an applied magnetic 
field, indicated a non-conservation of Andreev states^ 3 .. 
Further theoretical work would be necessary to clarify 
this issue. Experimentally, the technique presented in 
our present work, measuring the evolution of the An- 
dreev states under a supercurrent, may in fact provide 
a definitive way to distinguish between the specular and 
non-specular cases^ 4 ^. 

Finally, we note that ZBCP were also seen on {001} 
samples, showing the same suppression effect by currents 



applied predominantly along <100>. Figure [3| plots an 
example of such ZBCP at 4.2K. Such ZBCP have previ- 
ously been observed on c-axis cuprate thin films 4 ^, and 
could be explained by the presence of terrace edges on 
the film surface (see inset of Fig0J which would allow 
in-plane tunneling to occur on nominally c-axis sam- 
ples. The probability for in-plane tunneling could also 
be enhanced by the inherently larger matrix element for 
in-plane versus c-axis tunneling 46 . Consistent with this 
picture are instances where the ZBCP is seen together 
with a gap structure at the same tip location on {001} 
samples. The main panel in Fig0| shows such a spec- 
trum, with a ZBCP flanked by well-developed gap edges 
at w±14meV. The gap structures could be identified as 
the rf-wave averaged density of states associated with 
c-axis tunneling 6,45,46 . Since c-axis tunneling does not 
involve Andreev interference, the gap structures are ex- 
pected to depend only on the OP amplitude and be rela- 
tively insensitive to the OP phase fluctuations. What we 
observe in FigQ] is entirely consistent with this expecta- 
tion. Upon the application of a supercurrent, the ZBCP 
is suppressed as before but without significant variation 
in the gap structure. This is a remarkable observation, 
indicating that an applied supercurrent could weaken the 
local phase coherence while maintaining the pairing am- 
plitude. 

V. CONCLUSION 

In summary, using a scanning tunneling spectroscopy 
technique we have observed nanoscale suppression of low- 
energy Andreev states by a directly-applied supercurrent 
in YBa2Cu307_,5 thin-film strips at 4.2K. The distinct 
non-conservation of Andreev states indicates that this 
phenomenon cannot be easily attributed to extrinsic ef- 
fects such as Joule heating, vortex flow or excess quasi- 
particles. We have analyzed this suppression in terms of 
order-parameter phase fluctuations, using a model based 
on the af-wave BTK formalism. Qualitative agreement 
between our model and data suggests that the super- 
current may induce local phase fluctuations in the high- 
temperature superconducting cuprates. Further work is 
needed to elucidate this nanoscale phenomenon. 
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